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We study the magnetic field dependence of muonium-antimuonium conversion 
induced by neutral (pseudo) scalar bosons. Only the SS operator contributes to 
the conversion of polarized muonium, but it gets quenched by a magnetic field of 
strength 0.1 Gauss or stronger. Conversion induced by SS couplings for unpolarized 
muonium is independent of magnetic field. Magnetic fields of 0.1 Tesla or stronger 
starts to suppress conversion induced by PP interactions in the lowest Breit-Rabi 
level, but gets partially compensated by a rise in conversion probability in the other 
unpolarized level. The effects of [S^P){S^P) and {S^P){S ±P) operators behave 
in the same way as {V ^ A){y ^ A) and {V ^ A){V ^ A) operators, respectively. 



The spontaneous conversion of muonium (hydrogen-like atom M = fi'^e~) into antimuo- 
nium (M) would violate the separate additive muon and electron numbers, but would remain 
consistent with multiplicative muon or electron number conservation. Defining the effective 
coupling via the interaction 

'Hmm = ~^^^'y^(^ ~ 75)e/i7^(l - 75)e + h.c, (1) 

the limit has just been improved by an order of magnitude, 

Gmm < 1-8 X 10-2 Gf, (2) 

compared with the previous bound of 0.16 Gp 0, where Gp is the Fermi constant. The 
ultimate aim is to reach the sensitivity level of 10^^ Gp. 

The {V — A)(y — A) interaction of eq. (1) gained further theoretical footing when Halprin 
I^J^ pointed out that in left-right symmetric models with Higgs triplets, doubly charged 
scalars can mediate M-M transitions at tree level. The effective interaction is of {V ± 
y4)(\^ ± A) form after Fierz rearrangement. The possibility of {V — A){V + A) interactions, 
induced by dilepton gauge bosons P, was discussed by Fujii et al. [j^. Interestingly, muonium 
conversion is more pronounced in the singlet channel, in contrast to the {V — A){y — A) 
case where the conversion matrix element is of equal strength for both singlet and triplet 
muonium [jl|. 

Recently, Horikawa and Sasaki pointed out further that (V — A){y + A) interactions 
are considerably less senstive to magnetic fields compared to the (V — A)(y — A) case. 
Since actual experiments are conducted in the presence of magnetic fields this has 

important implications for the interpretation of the bound of eq. (2) for {V — A){y + A) 



interactions. In a recent paper |T0| (see also |TT| , p!2| ), we explored neutral (pseudo) scalar 
induced M-M oscillations. Since (5 =F -P) (^ ± P) and (V =F ^) (^ ± A) operators are related 
by Fierz transform, we wish to study the magnetic field dependence of neutral scalar induced 



muonmm conversion. 



Exotic neutral (pseudo) scalar bosons H and A could have couplings |T0 



— Cy = fie H + i-^ fi'y^e A + h.c, (3) 
V 2 v 2 



while a discrete symmetry such as multiphcative electron number Pg [1T3| forbids processes 
odd in number of electrons (plus positrons) like /x — 67 and /i eee. The resulting effective 
Hamiltonian responsible for M-M conversion is, 



^MA? = ~ 7r\ P'l5e /iTse 



where 5 = /le and P = fi-f^e. Note that 5^ - P^ and + P^ contain {S =F P)(S' ± P) and 
{S =F -P)('S' =F -P) interactions, respectively. In the "f/(l) limit" of fn = /a and m/^- = 111^ 
the subleading (5 =F P)(5' =F -P) terms are completely absent. 

The matrix elements of eq. (4) and the accompanying phenomenology of eq. (3) have 



been discussed in ref. [0. To explore magnetic field dependence of muonium conversion 



probabilities, consider the Hamiltonian for IS* muonium, 

n = Ho + a Se ■ Sp - ij,^ ■ B - ij,p^ ■ B + Hmm, (5) 

where Tio gives the IS* energy Eo = — a;^m/2, with reduced mass l/m = l/rrie + ^/fn^, 
a = 1.846 X 10~^ eV is the IS' muonium hyperfine splitting, and /Xg = —g^fiBSe, M/i = 
+9^il^B^ s^' where = g^, = 2 and /i^ = e/(2me) = 5.788 x 10~^ eV/Gauss is the Bohr 
magneton. Introducing the dimensionless parameters, 

x.Y^tsELi^n^A, (6) 

a \ rn^ ) 

we see that |y| is just 1% smaller than |X|. Ignoring Hmm for the moment, the four 



muonium Breit-Rabi energy levels IT^ 



Em{1, 0)=Eo + ^i^-^ + VTTX^ 

Em(0, 0)=Eo + |(-i-v^rTX^), (7) 



which correspond to the eigenstates 



|M;1,±1)= |M;n), |M;U), 
|M;1, 0)= c|M;n) + s|M;iT), 
|M;0, 0) = -s|M;Ti) + c|M;iT), 



(8) 



where the magnetic field dependent "rotation" is 
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(9) 



We have labeled the Breit-Rabi energy levels with the weak field basis \M] F^mp) , i.e. 
the corresponding zero B field {X,Y — > and s, c — > 1/^2) hyperfine states. The 
"uncoupled" basis of \M\'ms^,msp) corresponds to the strong field hmit of X, Y ^ 1 and 
s — ^ 0. In this limit, |M; 1, 0) — > |M; ti) and |M; 0, 0) — > |M; it), and the electron and 
muon hyperfine spin-spin coupling is overwhelmed by the Zeeman effect. 

For antimuonium, again ignoring the effect of H-mmi retaining the spin labels and with 
uncoupled basis |M;ms^,m5^), one simply fiips X — > —X, Y — > —Y and hence interchange 
s <-> c in eqs. (8) and (9). The notable changes are 



£;m(i,ti) = £;m(i,±i), 

since for given spin, the antiparticle magnetic moments have fiipped sign, and, 

|M;1,0)= s|M;n) + c|M;iT), 
|M;0,0) = -c|M;n>+s|M;iT>- 

We have the following energy differences between M and M eigenstates, 
EM{l,±l)-Ej^{l,±l)^±aY, 

Em{1, 0)-Em{1, 0)= Em{0, 0)-Em{0, 0) =0, 

Em{1, 0)-Em{0, 0)^-{Em{0, 0)-Em{1, 0))^aVlTX^. 



(10) 



(11) 



(12) 



The effect of Hmm can be treated as a perturbation. Define generically 



(M|Hmm|M) = 5/2 (13) 

(for simplicity, we take 5 to be real) between any two Breit-Rabi M and M energy eigen- 
states, it was shown by Feinberg and Weinberg |T]] that the time integrated probability for 
an initial muonium state to decay as antimuonium is 

where A = Em — is the energy difference, and A = 2.996 x 10"^° eV is the muon decay 
rate. The physics is clear: muonium oscillation has to compete with muon decay and the 
damping from oscillations between two states that are very disparate in energy. The total 
transition probability is 

Pt{M)= \cF,n.A' P^^''^^\M) (15) 

F,mp 

where |ci?^mf.p are the populations in muonium states of eq. (8), and p(-^''"^)(M) are the 
probabilities for an initial {F,mp) muonium state to decay as antimuonium. 

In principle P^^'^\M) and P^^''^\M) also contain the probabilities of initial (1, 0) or (0, 0) 
muonium states to decay as (0, 0) or (1, 0) antimuonium, respectively, since (1, 0) and (0, 0) 
are mixtures of unpolarized states. To show that these are vanishingly small, it is useful to 
notice the hierarchy 

5 < A < a, (16) 

the first of which follows from eq. (2) for any T^mm model. Combining eqs. (12) and (14) we 
see that (1, 0) (0, 0) and (0, 0) (1, 0) transitions are extremely suppressed by hyperfine 
splitting even for B = 0, and need not be considered. Similarly, although (1,±1) (1,±1) 
transitions contribute in i? = limit, they become rapidly suppressed even for rather weak 
magnetic fields because of the mismatch between Zeeman energy levels for M vs. M 
polarized states with same mp. 

In this work we will discuss only the relative magnetic field dependence of M-M conver- 
sion probabilities. Differences in coupling strength for various effective operators at zero B 



field can be found in refs. ||5|-[7|,[T0| . Hence, we normalize effective interactions to the conver- 
sion probability due to eq. (1) at zero magnetic field, and in particular for equally populated 
(1/4 each) Breit-Rabi states (the latter condition would be removed at the end). Thus, the 
zero field total transition probability is 

P^M; = 0) ^ ^ = 2.56 X 10-^ (^) ' , (17) 

which defines the parameter 5. 

For the {V — A){y — A) case one basically multiplies the r.h.s. (right hand side) of eq. 
(13) by 5ra,_^msMs^,ms-^^ ^euce 

(M;1,±1|7^mm|M;1,±1) = ^, 
(M; 1, 0|Hmm|M; 1, 0) = (M; 0, mMuW-, 0, 0) = ^^^TX^ ' ^^^^ 



and all other matrix elements vanish. One therefore finds the result []I5 



p(i,±i)(]Vf) 



2(^2 + a2F2 + A2) ' 

P(l' 0).^)^ p(0, 0).^N _ /^gX 

^ ' ^ ' 2(52 + A2(l + X2)) ^ ' 

The magnetic field dependence for Pt{M), as well as the separate probabilities P^^'^^^^M), 
P(^'°)(M) and p(°'°)(M) of eq. (19), are plotted as "+" symbols in Figs. 1-4, respectively. 
The behavior is readily understood. Because of the aY energy splitting, the suppression of 
(1, ±1) modes sets in with B field of just a few cG, and they become quenched for 0.1 G 
or higher. The mi;' = "unpolarized" modes are oblivious to the magnetic field until X 
becomes appreciable, i.e. for B ~ a/2fiB ~ IkG, and get quenched by fields of 1 Tesla or 
higher. The scale difference for Fig. 4 would be discussed shortly. The suppression in P(M) 
has been taken into account in the experimental limit of eq. (2) for the interaction of eq. 
(!)• 

We have checked and confirmed the result for the {V — A){V + A) case 0, 

52 



p{i,±i)(M) 



6(52 + a2F2 + A2) 
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(20) 



The results are also plotted in Figs. 1-4 as "x" symbols. Conversion in (0,0) mode is the 
most prominent 0, but gets suppressed by up to a factor of 4/9 when the magnetic field 
goes beyond ~ IkG. The (1, ±1) modes are quenched by magnetic fields of O.IG or higher, 
just like in the {V — A)(y — A) case, but P^^'^\M) actually grows with B field around IkG, 
and partially compensates for the drop in p(°'0). 

We now state the results for the (pseudo) scalar induced interaction case. Details would 
be given elsewhere [|T6|. For purely scalar interactions (/^ = 0), we find 

pa'±i)(M) = ^- , 

o)(^)^p(o, °)(M) = ^^^i^. (21) 

Thus, aside from the familiar quenching of the (1, ±1) states, the (1, 0) and (0, 0) states are 
completely insensitive to magnetic fields [|12|. For purely pseudoscalar interactions {fn = 0), 



we find 

p(i'±i)(M) = 0, 



2 2 



(22) 



In this case, muonium conversion occurs solely in the (0, 0) mode for zero magnetic field 
T0| . Conversion in the (1,0) mode starts to grow from zero for field strengths beyond ~ 
IkG, and partially compensates for the drop, by a factor of 4, in transition probability in the 
(0, 0) mode. We plot the results again in Figs. 1-4, with sohd and dashed lines representing 
scalar and pseudoscalar case, respectively. 

The results for {S — P){S — P) and {S — P){S + P) operators can be similarly obtained. 
With the same normalization conditions as described above, the results are also given in 
Figs. 1-4 for sake of comparison, with open circles representing the {S — P){S + P)case 
and open boxes representing the {S — P){S — P) case. Note that according to eq. (4), the 



{S—P){S—P) operators should be subdominant compared to {S—P){S+P) operators. Pure 
{S =F P)(S' ± P) operators corresponds to complex neutral scalars [|1^], where the sneutrino 
Ur in SUSY models with i?-parity breaking [|T^ as a special case. It is evident from Figs. 
1-4 that the combinations {S — P){S — P) and {S — P){S -\- P) behave in the same way as 
(y — A){V — A) and {V — A){V + A), respectively. In the latter case, the two operators are 
related to each other by a Fierz transform. For the former case, although the operators can 
not be related to each other by a Fierz transform, the matrix elements are always in same 
proportion, which comes as a consequence of the nonrelativistic limit. 

Turning to discussions, we note that the assumption of equally populated Breit-Rabi 
levels is not a valid one, since this is determined by the muonium formation process and 
the magnetic field strength. However, as a consequence of this assumption, the results in 
eqs. (19-22) and hence Figs. 2-4 all have an artificial factor of 4. To illustrate the effect 
of differently populated Breit-Rabi levels, normalizing again to the {V — A){V — A) case 
at zero B field, we take muonium states to be populated as [^] 32%, 35%, 18% and 15%, 
respectively, for (F, mp) = (0, 0), (1,+1), (1, 0), (1, -1), and plot the resuhs for Pt{M) 
(i.e. analogous to Fig. 1) in Fig. 5. The (V — A)(y — A) and {S — P){S — P) results are 
unchanged, since the conversion matrix elements are the same for all modes. The purely 
scalar (SS) case is also unchanged, since |cioP + |cooP = |ci+ip + |ci_ip = 50% is the 
same as the equally populated case. The (V — A)(y + A), {S — P){S + P) and PP cases 
are somewhat modified from the equally populated case, but the difference for IkG field is 
rather slight. 

Let us summarize our findings. Purely scalar (SS) induced M-M transitions in polarized 
(i.e. {F^mp) = (1,±1)) modes are quenched by magnetic fields of 0.1 Gauss or higher, 
but conversion in unpolarized states (i.e. (1,0) and (0,0)) are independent of magnetic 
field strength. Purely pseudoscalar (PP) interactions do not induce conversion in (1,±1) 
states, but exhibit compensating effects in the (1, 0) and (0, 0) channels, similar to the 
{V-A){V + A) case. Interactions of {StP){S±P) and {VtA){V±A) form have the same 
magnetic field dependence since they are related by Fierz transform, while {S =F P)(S' =F P) 

8 



and {V ^ A){V ^ A) interactions have the same field dependence because of proportional 
conversion matrix elements. 
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FIGURES 

FIG. 1. Magnetic field dependence of total muonium conversion probability Pt{M) assuming 
|ci,ip = |ci,oP = |ci,_ip = |co,oP = 1/4, and normalized to conversion strength of {V — A){V — A) 
interaction at zero magnetic field. Solid and dashed lines stand for SS and PP operators, re- 
spectively, while o, □, + and x stand for {S - P){S + P), {S - P){S - P), {V - A){V + A) and 

iy — A){V — A) cases, respectively. 

FIG. 2. P^^'^^\M) vs. magnetic field with same assumption as Fig. 1. 

FIG. 3. P^^''^\M) vs. magnetic field with same assumption as Fig. 1. 

FIG. 4. P^^'^'^iM) vs. magnetic field with same assumption as Fig. 1„ 
FIG. 5. The same as^ Fig. 1 except |ci,iP = 0.35, |ci,or = 0-18, |ci -ip = 0.15 and 

|cooP = 0.32. 
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